Magnetic Properties of the Second Mott Lobe in Pairing Hamiltonians 
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We explore the Mott insulating state of single-band bosonic pairing Hamiltonians using analytical 
approaches and large scale density matrix renormalization group calculations. We focus on the sec- 
ond Mott lobe which exhibits a magnetic quantum phase transition in the Ising universality class. 
We use this feature to discuss the behavior of a range of physical observables within the framework 
of the ID quantum Ising model and the strongly anisotropic Heisenberg model. This includes the 
properties of local expectation values and correlation functions both at and away from criticality. 
Depending on the microscopic interactions it is possible to achieve either antiferromagnetic or ferro- 
magnetic exchange interactions and we highlight the possibility of observing the Es mass spectrum 
for the critical Ising model in a longitudinal magnetic field. 

PACS numbers: 37.10. Jk, 75.10.Jm, 75.10.Pq 
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I. INTRODUCTION 

The observation of Bose-Einstein condensation (BEC) 
in dilute alkali gases [l|, Q has led to a wealth of activity 
linking ultra cold atoms research and condensed mat- 
ter physics. The precise control over atomic interactions, 
and the use of optical lattices, offers valuable insights into 
the effects of strong correlations Q. This is exemplified 
by experiments on the Bose-Hubbard model which reveal 
the quantum phase transition from a superfluid (SF) to 
a Mott insulator (MI) as the depth of the optical lattice 
is increased [J, [f|. Motivated by these advances, recent 
attention has been directed towards a variety of multi- 
component systems, including spinor condensates [6T4l7| 
and atomic mixtures [l8rl23j . The possibility of novel be- 
havior is greatly enhanced in the presence of these addi- 
tional "spin" degrees of freedom, and routes to quantum 
magnetism have been proposed by exploiting internal hy- 
perfine states [24[ and using different atomic species pH - 
|28[ . The decoupling of the electronic and nuclear spins in 
alkaline earth atoms has also been suggested as a way to 
realize quantum spin liquids and exotic magnetism based 
on the SU(N) groups [2§-[3l|]. 

In recent work some of the present authors suggested 
using atom-molecule mixtures as a route to the paradig- 
matic quantum Ising model [lU, ■ Motivated by stud- 
ies of the BEC-BCS transition for bosons, both in the 
continuum limit [34H36! ] and on the lattice (37l - l40| . we 
investigated the rich phase diagram of bosons interacting 
via Feshbach resonant pairing interactions in an optical 
lattice. Combining exact diagonalization (ED) on small 
systems with analytical predictions based on the strong 
coupling expansion, we provided evidence for an Ising 
quantum phase transition in the second Mott lobe. In 
contrast to previous numerical studies which advocated 
the presence of super-Mott behavior [39|, H(| , the Ising 
spectral gap indicates the absence of low l ying sup erfluid 
excitations deep within the Mott phase [32|, [33|- This 



conclusion gained further support in a recent comment 
PH which shows the presence of exponential decay in the 
connected correlation functions. A detailed discussion of 
the superfluid properties [34! - 36] in one dimension (ID) 
was also provided in Ref. |42| using large scale density 
matrix renomalization group (DMRG) [43[ and field the- 
ory techniques. 

In view of the broad interest in these systems, and the 
technical issues surrounding Refs. [H,[4(|, we discuss the 
properties of the second Mott lobe in detail using analyt- 
ical arguments and DMRG. To aid the comparison with 
these previous works we focus primarily on the homonu- 
clear case with a single species of bosonic atom. We also 
begin by restricting the local Hilbert space for the atomic 
and molecular occupations In this reduced setting 
we determine both the locus of the antiferromagnetic 
Ising transition and the onset of superfluidity. We also in- 
vestigate the atomic and molecular correlation functions 
within the MI and compare directly to the predictions of 
the quantum Ising model. As advocated in Refs. [H, HH 
this provides a simple and intuitive framework in which 
to discuss the absence of super-Mott behavior (39l - [4l| . 

Going beyond the restricted Hilbert space description 
we show that the same ideas apply. Interestingly, by 
tuning the microscopic interactions it is also possible to 
change the sign of the Ising exchange interaction from an- 
tiferromagnetic (AFM) to ferromagnetic (FM). In general 
these Isin g H amiltonians also contain an effective mag- 
netic field [321.1331] and in ID this will act as a confinement 
potential for the zero field excitations of the FM chain 
HI, SI] • This suggests the possibility of observing the 
non-trivial Eg mass spectrum of "meson" bound states 
for the critical FM Ising model in a longitudinal field 
[iril - |48j . As an extension of these results we also provide 
the magnetic Hamiltonian for Bose-Fermi mixtures. 

The layout of this paper is as follows. In Sees. Inland 
IIIII we provide an introduction to the bosonic Feshbach 
Hamiltonian and discuss the mapping to the quantum 
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Ising model. In Sec. II VI we present a cross section of the 
phase diagram obtained by DMRG which displays both 
the Ising quantum phase transition and the onset of su- 
perfluidity. We confirm the Ising behavior within the 
Mott phase using results for the excitation gap and the 
entanglement entropy. In Sec. [V]we investigate the role 
of higher order terms in the strong coupling expansion. 
We discuss their impact on the non-universal properties 
of the phase diagram such as the curvature of the Ising 
phase boundary. In Sec. lVII we examine the local expecta- 
tion values within the Mott phase at and away from crit- 
icality. We move on to correlation functions in Sees. IVIII 
and I Villi and discuss their relation to the Ising model. In 
Sees. IIXI and Ixl we turn our attention to softcore bosons 
and demonstrate the existence of both antiferromagnetic 
and ferromagnetic Ising transitions. We comment on the 
closely related fcrmionic problem in Appendix [B] Wc 
conclude in Sec. IXII and provide perspectives for further 
research. 



II. MODEL 
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We consider the Hamiltonian j37H40j 

H = '^e a n ia - '^ S ^t a {b\ a b ja +H.c.) 

ioL (ij) a 

+ ^2 ~Y~ '■ n i a n la ' ■ +H F , (1) 

iaa' 

describing bosons, bi a hopping on a lattice with sites i 
where a = a,m labels atoms and molecules and rn a = 
b\ a bi a . Here, e a are on-site potentials, t a are hopping 
parameters, (ij) denotes summation over nearest neigh- 
bor bonds and U aa i are interactions. We use normal or- 
dering symbols to indicate : rii a ni a > : = rii a (ri{ a — 1) for 
like species and : ni a rii a > : = ni a rii a r for distinct species. 
For simplicity we assume that molecules are formed by 
s-wave Feshbach resonant interactions 



g y^X m l a i a i + H.C.), 



(2) 



where m, = bi m and a, = &ja! for recent work on 
the p-wave problem see Ref. [49(. An important fea- 
ture of the Feshbach interaction ^ is that atoms and 
molecules are not separately conserved. However, the to- 
tal, At = J2i( n ia + 2n ?I „) is preserved. One may there- 
fore work in the canonical ensemble with px = Nt/L 
held fixed, where L is the number of lattice sites. In or- 
der to make contact with the previous literature (39M4l| 
we consider only a single-band description in E g . (HI) ; for 
a discussion of higher band effects see Refs. [5XM32I ]. As 
reported elsewhere the pairing Hamiltonian ([1]) has a rich 
phase diagram exhibiting both MI and SF phases [32| — 
liol ■ Most notably, the system displays a discrete Z 2 
symmetry breaking transition |34| - l36f between a paired 
molecular condensate (MC) and an atomic plus molecu- 
lar condensate (AC+MC) phase [42j; for closely related 
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FIG. 1. (color online). Depiction of the second Mott lobe 
of the Hamiltonian {TJ with px = 2 showing either a pair of 
atoms or a molecule on each site. These may be regarded as an 
effective spin-down (JJ.) and spin-up (ff) respectively. Second 
order virtual hopping processes lead to J zz interactions and an 
effective Ising model. XY exchange, J xy , occurs at third order 
and involves interchanging a pair of atoms and a molecule. 
The Feshbach term @ converts a pair of atoms (4) into a 
molecule (ff) and therefore acts like a transverse field. 



transitions in other models see also Refs. [53M57I |. Here, 
our main focus will be on the MI phase. In particular, we 
shed further light on the magnetic characteristics of the 
second Mott lobe [HI, HU . We also discuss the connec- 
tion to Ref. and earlier numerical studies H^, |4(j ■ In 
order to facilitate our numerical simulations we consider 
the restricted Hamiltonian 



Hp = THV, 



(3) 



where the projection operator V projects on to the re- 
stricted local Hilbert space with a maximum of r a atoms 
and r m molecules per site. We begin our discussion in 
Scc.|TTT]by setting r a = 2 and r m = 1 as used in Ref. [4fj| . 
As we will see, the essential characteristic features of the 
phase diagram are readily gleaned from this limiting case. 
We move on to the more general problem with canoni- 
cal softcore bosons in Sees. IIXI and IXl Throughout this 
manuscript we use the value of U am to set the overall 
energy scale; in Figs. [2] [3l |4] and Figs. El [3 [8] we set 
U am = 1 , and in the remaining figures we set U am = 4 in 
order to descend deeper in to the Mott phase. 



III. SECOND MOTT LOBE 

A convenient way to describe the second Mott lobe 
with px = 2 is via an effective spin model derived within 
the strong coupling expansion; see Fig. [TJ Introducing 
effective spins | if) = |2;0)/\/2 and | ft) = |0; 1) in the 
occupation basis \n a ;n m ) one obtains the effective spin- 
1/2 quantum Ising model [321 [33j| 

H ~ J zz Y^S!Sf +1 +hJ2S?+Tj2sf+C+0(t 3 ), (4) 

i i i 

where we work to second order in the hopping parame- 
ters. The spin operators are given by 

5+ = mW s _ 



aVm _ (n m - np/2) 

and S ± = S x ± iS y . The factors of \[2 arise from 
the action of the Bosc operators on the basis states. 
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Equivalcntly, since n a + 2n m = 2, one may also write 
S z = (1 - n„)/2 or S z = (2n m - l)/2. As indicated 
in Fig. [IJ the first term in Eq. ((4]) arises from the vir- 
tual hopping processes of atoms and molecules on to 
neighboring sites and corresponds to an effective mag- 
netic exchange interaction, J zz . The second term in 
Eq. (0} reflects the energetic detuning between atoms 
and molecules and corresponds to an effective longitu- 
dinal magnetic field, h. The third term in Eq. (|4]) corre- 
sponds directly to the Feshbach term in Eq. $2$. Indeed, 
it is readily seen from the definition ([5]) that S + converts 
two atoms (JJ.) into a molecule (ff) and therefore acts like 
a spin raising operator. It follows that the Feshbach term 
(l2l) acts like a transverse field in the x-direction with 



T = 2gV2. 



(6) 



The overall structure of the magnetic Hamiltonian (|4]) 
is generic to the second Mott lobe with px = 2. How- 
ever, the coefficients depend on the specific Hilbert space 
restriction. For the Hamiltonian ([3]) with restriction pa- 
rameters r a = 2 and r m = 1 one obtains [32l . [Hf 



4/ 



U a 



u a 



U a 



and 



h = e m - 2e a - U aa . 
The constant offset is given by 
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(7) 



(8) 



(9) 



where z is the cubic lattice coordination and z = 2 in ID. 
In general there will also be additional contributions to 
these coefficients arising from higher order terms in the 
strong coupling expansion. As we will discuss in more 
detail below similar results also hold for canonical soft- 
core atoms and molecules but with modified coefficients. 
This is due to the presence of additional intermediate 
states that are explored in the virtual hopping processes. 
Nonetheless, this truncation of the Hilbert space is useful 
for initial numerical simulations, and is also used in Rcf. 
[40j . In addition the principal features of the magnetic 
description are more readily exposed. We will return to 
canonical softcore bosons in Sees. HXl andlXl 



track the Ising quantum phase transition throughout the 
second Mott lobe, or to accurately delimit the onset of 
superfluidity We address this situation for the homonu- 
clear case by using DMRG on large systems. We keep up 
to 3000 density matrix states so that the discarded weight 
is less than 10~ 10 . Following Rcf. [4(J we implement the 
model ^ as a two-leg ladder system, where the atoms 
and molecules reside on opposite legs of the ladder. In 
this representation the Feshbach term ([2]) corresponds to 
hopping along the rungs, and extreme care must be taken 
for small values of g. In Fig. Owe present a cross sec- 
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FIG. 2. (color online). Phase diagram of the ID Hamiltonian 
© with p T = 2 obtained by DMRG with up to L = 256 
sites. We restrict the local Hilbert space to a maximum of one 
molecule and two atoms per site with r m = 1 and r a = 2, and 
set e a = e m = 0, U aa = 0, U am = 1, corresponding to h = 0. 
We set t a — 2t m = t and vary t to obtain different values 
of Jzz = 17£ 2 /4 corresponding to antiferromagnetic (AFM) 
exchange. This cross section shows a Z2 ordered AFM MI, 
a Z 2 disordered AFM MI, and an AC+MC SF. The filled 
squares are obtained from the vanishing of the first excitation 
gap, Ai = E\ — Eo corresponding to an Ising quantum phase 
transition within the MI; the solid line is a spline fit which is 
extrapolated down to T — 0. The dashed line corresponds to 
the strong coupling result, T c = J zz /2 which follows from the 
Hamiltonian (Q with h — 0. The dotted line is obtained by 
ED of Hamiltonian (T3]) including the subleading XY exchange 
terms, J xy = —0.46 Jzl 2 which arise at higher order in the 
strong coupling expansion. The diamonds are obtained from 
the vanishing of the two-particle gap, E2 g = fi+ — fi- where 
H± = ±[E {L,N T ± 2) - E a (L,N T )} and up to L = 128, 
indicating the onset of an AC+MC superfluid. 



IV. PHASE DIAGRAM AND THE 
ANTIFERROMAGNETIC ISING TRANSITION 

In Refs. HH, [33| we provided a variety of evidence for 
an Ising quantum phase transition in the closely related 
heteronuclear generalization of Hamiltonian ([1}. How- 
ever, due to the small system sizes accessible by exact 
diagonalization, a complete elucidation of the phase dia- 
gram was not possible. In addition correlation functions 
were out of reach. In particular, it was not possible to 



tion of the phase diagram for the Hamiltonian with 
r a = 2 and r m = 1. In this approach the Hilbert space 
restriction parameters r a and the interactions U aa > are 
treated independently. For ease of exposition we begin 
by setting e a = e m = U aa = and U am = 1, where our 
choice of parameters is motivated by the simplest case 
with zero magnetic field as given by Eq. ||5J). In addi- 
tion, for small hopping parameters t a , and large U am , 
one may set U aa = in Eq. (J7J) without loss of gener- 
ality or conflict with the strong coupling expansion. We 
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will incorporate the effects of finite U aa in our subsequent 
discussion. 

The phase diagram in Fig. [2] contains three distinct 
phases. A Z2 disordered MI with vanishing staggered 
magnetization J2i(~ l) l {Sf)/L = 0, a Z 2 ordered MI 
with a finite staggered magnetization and long range an- 
tifcrromagnctic correlations, and an AC+MC supcrfluid 
with both atomic and molecular power law superfluidity. 
The additional MC phase with only molecular superflu- 
idity is absent in this cross section of the phase diagram 
due to our choice of parameters; for more details of the 
superfluid phases see Ref. (42[. In this manuscript our 
main focus is on the MI phase. In particular, we see that 
the magnetic phase boundary bends over quite consid- 
erably due to higher order terms in the strong coupling 
expansion. Nonetheless, the quantum phase transition 
remains in the Ising universality class. For example, the 
Ising character of this transition is supported by Fig. [3] 
which shows the characteristic linear variation of the ex- 
citation gap HU 

A = |r-r c |, (10) 

on passing through the transition in the scaling regime. 
The Ising character is also confirmed by DMRG results 
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FIG. 4. (a) DMRG results for the entanglement entropy, 
Sl(1), at the Ising transition within the MI for a subsystem of 
length I in an periodic chain with L = 64. We use the same 
parameters as in Fig. [2] and set t = 0.17 corresponding to 
J zz « 0.12 and T c « 0.028. The solid line is a fit to Eq. pi]). 
The extracted value of the central charge c = 0.50(2) is consis- 
tent with a magnetic transition in the Ising universality class, 
(b) The entanglement entropy difference AS(L) for t — 0.17 
shows a peak at the Ising transition. The peak height corre- 
sponds to c = 1/2 as indicated by the dashed lines. 



where s\ is a non-universal constant and c is the central 
charge. As shown in Fig.[U[a) the results are in excellent 
agreement with an Ising quantum phase transition with 
c = 1/2. In particular, this continues to hold further out 
in the Mott phase where the strong coupling analysis no 
longer strictly applies. Further evidence for this may be 
seen from the entanglement entropy difference [42l . |63| 



AS(L) = S L (L/2)-S L/2 (L/4) 



(12) 



FIG. 3. DMRG results for the excitation gaps, Ai = E\ — 
Eo (circles) and A2 = E2 — Eo (squares), for the bosonic 
Hamiltonian ((3| with open boundaries and up to L — 256 
sites. We pass through the Ising quantum phase transition 
shown in Fig. [2]with t — 0.05 and J zz w 0.01. The transition 
occurs at F c ~ 4.6 x 10 -3 which is slightly below the strong 
coupling result, F c = J zz /2 « 5.3 x 10~ 3 . This reflects the 
curvature of the Ising transition shown in Fig. The solid 
lines indicate the linear gap, A = |T — T c ], corresponding to 
the Ising critical exponent, v = 1. 

for the entanglement entropy as shown in Fig. [4] For a 
block of length / in a system of length L, the von Neu- 
mann entropy is given by Sl{1) = — Tr; (pjlnp;), where 
pi = Tr£_;(p) is the reduced density matrix. In a critical 
system with periodic boundaries one obtains [59l - [62| 



Sl(1) = -In 



"2L , 




— sin 




7T ' 





(11) 



which exhibits a peak on transiting through the mag- 
netic transition as shown in Fig. 0|b). At criticality 



AS = § In 2 + . . . [59H62J and the peak height is in good 
agreement with c = 1/2. It is interesting to note that 
a naive spline fit to the currently available DMRG data 
for the Ising phase boundary shown in Fig. [5] terminates 
within the ML This tentatively suggests that a magnetic 
transition may exist even in the absence of the Feshbach 
term ([2]). Unfortunately, it is difficult to gain a quan- 
titative handle on the character of this feature, which 
occurs for intermediate hopping strengths and vanish- 
ing Feshbach coupling, using our current implementation 
of DMRG. Moreover, this T = limit differs somewhat 
from other investigations of the two-component Bose- 
Hubbard model [ij, [65| , since we keep the total density, 
px = ^2ii n ia + 2rii m )/L = 2 held fixed, with two atoms 
and one molecule per site. It would be instructive to 
explore this limit in future work. 



V. HIGHER ORDER CONTRIBUTIONS IN 
THE STRONG COUPLING EXPANSION 

In the discussion above we have presented a variety of 
large scale DMRG results in favor of an Ising quantum 
phase transition taking place within the second Mott lobe 
[32l I33I ] . In particular we have argued that the Ising char- 
acter of the quantum phase transition persists through- 
out an extended region of the phase diagram as shown in 
Fig. [2] Nonetheless, it is important to bear in mind that 
the explicit Ising Hamiltonian given in Eq. Q is derived 
by means of the strong coupling, t/U expansion. It is 
therefore tailored towards a quantitative description of 
the bosonic Hamiltonian ([1} deep within the Mott lobe, 
as supported by our DMRG results. At larger values of 
the hopping parameters departures from the strong cou- 
pling result (|4]) are to be expected. This is evident from 
the deviation of the phase boundary from the asymp- 
totic result, r c = J zz /2 as shown in Fig. [5] Carrying out 
the strong coupling expansion to third order in the hop- 
ping parameters one must supplement the Hamiltonian 
(HI with the additional XY exchange terms: 

AH= J^j2(StSr +1 + S rSt +1 ), (13) 



where S ± = S x ± iS y , and 



J x 



4t z t 



1 



1 



Uam U aa V U am U am U aa 



(14) 



One thus obtains a strongly anisotropic XXZ Heiscn- 
bcrg Hamiltonian in both a longitudinal and transverse 
field. For the parameters chosen in Fig. [21 J zz = \7t 2 /A 
and J xy = — 4t 3 . This yields the XXZ anisotropy pa- 
rameter, S = J xy /J zz = — 16t/17 corresponding to 
6 ps -0.46^717. Within the MI shown in Fig. [2] this 
gives —0.27 < 5 < 0. The system remains in the strongly 
anisotropic regime and Ising criticality is expected to per- 
sist throughout this cross section of the MI. Nonetheless, 
the XY contributions (|13|) modify the location of the Ising 
quantum phase transition as shown in Fig. [2j As we will 
discuss in Sec. IVIII such terms also influence the non- 
universal amplitudes in the correlation functions whilst 
preserving the universal critical exponents. 
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c 

FIG. 5. (color online). Local expectation values deep within 
the Mott phase obtained by DMRG for the bosonic Hamilto- 
nian Q with periodic boundaries and restriction parameters 
r m = 1 and r a = 2. In order to demonstrate the broader va- 
lidity of our results, whilst suppressing higher order terms in 
the strong coupling expansion, we choose different parameters 
to those used in Fig. [2] We set e a = and e m = Uaa = 3.8, 
corresponding to h = and take U am = 4 and t — 0.005. 
(a) Staggered magnetization |(Sf )|. The theoretical result for 
the transverse field Ising model is given by Eq. (|15[) and is 
indicated by the solid line. Inset: The application of a small 
staggered field h s t yields the same results as in panel (a) in the 
limit h s t — 0. (b) Transverse magnetization (Sf). The solid 
line shows the Ising behavior given by Eq. (|16[) . The dashed 
line corresponds to {Sf ) = —1/ir which holds at criticality. 



VI. LOCAL EXPECTATION VALUES 

Having provided numerical evidence for an Ising quan- 
tum phase transition in the second Mott lobe of the 
Hamiltonian ([1} we turn our attention to the local ex- 
pectation values. The order parameter of the antiferro- 
magnetic transverse field Ising model @ is the staggered 
magnetization (— l) l (Sf). In the thermodynamic limit 
with h = and T < T c this is given by 6rj|67| 



where /3 = 1/8 is the Ising critical exponent; in the dis- 
ordered phase with T > T c the order parameter van- 
ishes. In order to verify this characteristic behavior in 
the bosonic Hamiltonian ([3]) we choose parameters deep 
within the Mott phase. As shown in Fig. 03a) our DMRG 
results arc in excellent agreement with Eq. (fT5)) . This 
confirms that any higher order corrections to the Ising 
description Q are small. In a similar fashion the trans- 
verse magnetization is given by (6(1 \0l\ 



Mm) = 



1 



(15) {S f) = - f g 



2r + J zz cos k 



2tt v /4r 2 + Jf z + ATJ ZZ cos k 



(16) 
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This dependence is confirmed in Fig. [5{b) both at and 
away from criticality. Note that at the critical point 
where T = T c the expectation value (Sf) = —1/ir is 
non-vanishing due to the transverse field. We have also 
checked the consistency of our DMRG results by applying 
a small staggered field h st to the bosonic Hamiltonian ([3]) , 
A£T 8t = -hst^i-iySf, where Sf is given by Eq. ©. 
In the limit h s t —J- this replicates the effect of sponta- 
neous symmetry breaking in our finite size simulations. 
The results converge to the same values as for h st = 0; 
see inset of Fig. E^a). 

Having established the validity of the explicit Ising 
Hamiltonian ([4} deep within the Mott phase, it is in- 
structive to see how this leading order behavior is mod- 
ified as one moves out towards the MFSF boundary. In 
Fig. inja) we show the evolution of the local magnetiza- 
tion (Sf) with increasing hopping parameters within the 
Z2 disordered phase shown in Fig. [2] Deep within the 
MI, (Sf) = as one would expect for the disordered 
phase of the transverse field Ising model @ with h = 0. 
However, at larger values of the hopping parameters a fi- 
nite uniform magnetization develops as indicated by the 
plateaus in Fig.[6]Ja). This is due to the presence of higher 
order terms in the strong coupling expansion. In general 
such contributions are expected to induce corrections to 
the leading order coefficients, J zz and C, and the mag- 
netic field, h, given in Eqs. O - ©■ However, it is evi- 
dent from Fig.[6jb) that this uniform contribution to the 
magnetization remains significantly below the saturation 
value of (Sf) = —1/2, indicating that the leading order 
description Q is a useful starting point in the broader 
phase diagram. Indeed, the staggered magnetization re- 
mains zero in the bulk of the system, as may be seen 
from the absence of oscillations in Fig. [(^a). As such 
we remain in the disordered phase of an Ising antiferro- 
magnet, albeit in the presence of an increasing uniform 
effective magnetic field. More generally, one may also 
transit through the Z2 ordered region shown in Fig. [5] 
As indicated in Figs.^c) and (d) this results in the onset 
of antifcrromagnctic oscillations in the local magnetiza- 
tion. In Fig.^e) we plot the evolution of the correspond- 
ing uniform and staggered magnetizations. The region of 
support of the staggered component is in agreement with 
the Z 2 ordered phase inferred from the gap data in Fig. 
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VII. CORRELATION FUNCTIONS 

Having discussed the local expectation values we now 
consider correlation functions. In order to orient the dis- 
cussion we first recall the expected theoretical behavior 
at the critical point of the antiferromagnetic Ising model 
in a transverse field, where T c = J zz /2 and h = 0. In the 
absence of a longitudinal field the asymptotic longitudi- 
nal correlation function decays as a power law [66| 

(S?S? +n )~(-irBn-" + ... J (17) 



FIG. 6. (color online), (a) DMRG results for the local magne- 
tization (Sf) for the parameters used in Fig. [2] with h = and 
T = 0.04. We set t = 0.025, 0.075, 0.125, 0.325 from top to 
bottom with L — 32 and periodic boundaries. The plateaus 
correspond to the development of a uniform magnetization 
whilst the staggered magnetization remains zero, (b) Uni- 
form magnetization Mu = ^2rASf) /L. The non-zero value 
is attributed to higher order terms in the strong coupling ex- 
pansion. These may modify the leading order magnetic field 
given in Eq. JSJ. The dotted line shows the approximate lo- 
cation of the MI to AC+MC transition obtained from the gap 
data in Fig. [2] (c) We set Y — 0.02 and use the same values 
of the remaining parameters as in panel (a) . The oscillations 
correspond to the onset of the Z2 ordered phase in Fig. [2] (d) 
Antiferromagnetic oscillations in the Z2 ordered phase with 
t = 0.125 (circles), t — 0.155 (squares), t = 0.175 (diamonds) 
and t = 0.195 (triangles), (e) Uniform magnetization A^u 
(circles), and staggered magnetization Als = (Sf) /L 

(squares) for Y = 0.02. 
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where T] = 1/4 is the Ising pair correlation exponent, 
B = _4~ 3 4" 1 2 1 / 12 e 1/4 ~ 0.161 is a constant prefactor, 
and A ~ 1.2824 is the Glaisher-Kinkclin constant. In 
a similar fashion, the transverse correlators also exhibit 
power law behavior at criticality [66j 



(SfSf +n ) ~ (Sf ) 2 + (2^n)- 2 + 
where (Sf) = — 1/tt at the transition and 



(SfSf 



i-\-ni 



-(-l) n (S/4) n" 9 / 4 



(18) 



(19) 



As suggested by equation (fT5)l , due to the finite value of 
(5f) = —1/tt at criticality one must consider the con- 
nected correlation function (SfSf +n ) — (Sf)(Sf +n ) in or- 
der to see power law behavior. To establish this depen- 
dence in the bosonic Hamiltonian ^ we perform DMRG 
calculations with open boundaries and up to L = 512 
sites. As shown in Fig. [JJ at T = T c the data are con- 
sistent with power law behavior in the bulk of the sys- 
tem away from the sample boundaries. In addition the 
data show long range order in \(SfS? +n )\ for T < T c , 
and a finite correlation length for T > T c . Although this 
provides direct evidence for a quantum phase transition, 
a quantitative determination of the Ising critical expo- 
nents requires a more detailed finite size scaling analysis 
of the data. This is most readily achieved using peri- 
odic boundary conditions. We recall that in a finite size 
critical system with periodic boundary conditions, the 
two-point function of a primary field O may be obtained 
by conformal transformation of the strip geometry (68j : 



(0(r 1 )0(r 2 )) L =M 



(20) 



Here a is the critical exponent in the thermodynamic 
limit and M is the prefactor: (O(ri)0(r 2 ))oo = Nr~ a . It 
follows from equation (f2"U)) that the rescaled combination 



L a (0(n)0(r 2 )) L = f a (r/L) 



is a prescribed scaling function 

fa(x)=M 



sin(7ra;) 



(21) 



(22) 



of the reduced separation x = r/L. As shown in Figs. [5] 
and [51 the rescaled critical two-point functions 



F z {n/L)=LV%SlS* +n ) L \, 
F y {n/L)=I?l%StSt +n ) L \, 
F x (n/L) ee L 2 [(SfSf +n ) (Sf)(Sf +n )}j 



(23) 



all show striking data collapse over the entire system 
length, showing clear indications of Ising criticality. Deep 
in the Mott phase the resulting scaling functions are in 
excellent agreement with the theoretical results for the 
lattice Ising model (U) in finite size geometry, as indi- 
cated in Fig. [S] This includes both the universal critical 
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FIG. 7. Correlation functions within the second Mott lobe of 
the bosonic Hamiltonian © obtained by DMRG with L = 512 
and open boundaries. We use the same parameters as in 
Fig. [3] The data are consistent with an Ising quantum phase 
transition between a magnetic state with long range order 
and a disordered state with a finite correlation length. The 
solid lines correspond to the critical two-point functions in 
Eqs. HjJ, (|18[) and (|19|) . A quantitative demonstration of 
the Ising critical exponents is shown in Fig. [8] using periodic 
boundaries and finite size scaling. 



exponents, a and the non-universal amplitude prcfactors, 
Af taken from Eqs. ([TTj> . ([T8]) and ([19]). Further out in the 
Mott phase, the non-universal prefactors are influenced 
by higher order terms in the strong coupling expansion 
as discussed in Sec. El but the universal Ising exponents 
are robust to these perturbations; see Fig. [8] 

Having confirmed the presence of a line of critical Ising 
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FIG. 8. (color online). DMRG results for the rescaled bosonic 
correlation functions L 1/4, \{Si Sf +n ) l (circles), L 2 [(Sf Sf +n ) — 
(Sf)(Sf +n }] L (squares) and L 9/4 |<Sj'£iV n )i,| (triangles) with 
periodic boundaries at criticality for the parameter set of 
Fig. [7] with t — 0.05. The data collapse over the entire system 
length with the Ising critical exponents. The non-universal 
prefactors differ slightly from the theoretical predictions of 
the lattice Ising model @ as indicated by the solid lines. This 
is due to the presence of small additional XY contributions 
to the Ising description. By descending deeper into the Mott 
lobe one obtains a complete quantitative agreement including 
the non-universal prefactors as shown in Fig. [9] 
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FIG. 9. (color online). DMRG results for the rescaled bosonic 
correlation functions L 1 ^ 4 ] (S* Si +n )L (circles), L 2 [(Sf Sf +n ) — 
(Sf)(S? +n }] L (squares) and L 9 ^\(SV SV +n ) l\ (triangles) with 
periodic boundaries at criticality for the parameter set of 
Fig. [5] with t — 0.005. The data show clear scaling collapse 
over the entire system length. The theoretical results for the 
lattice Ising model (Q in finite size geometry are indicated 
by solid lines. The data are in excellent agreement with both 
the universal Ising critical exponents and the non-universal 
amplitudes in Eqs. (Tl7j, CEH1> and (fl9)) . 



correlations within the second Mott lobe of the bosonic 
Hamiltonian (J3j) we now consider the generic behavior in 
the MI phase. As may be seen in Fig. [7J the data reveal 



lxio 




1x10 



FIG. 10. Off-critical order parameter correlations |(S'fS'f + „)| 
obtained by DMRG on the ID system ([3jl with L = 48 and 
periodic boundaries. We use the same parameters as in Fig. [5] 
The solid lines correspond to the quantum Ising model (|4]) and 
are given by Eq. §2%$ for T < T c and Eq. (J25J) for F > T c . 
The agreement confirms the presence of long range order for 
r < T c and exponential decay for F > T c . 



a finite correlation length for the connected correlation 
functions on either side of the transition. This is consis- 
tent with the presence of an Ising spectral gap as shown 
in Fig. [3] In order to place this behavior on a quantitative 
footing we recall the principal results for the transverse 
field Ising model (j4]) with h = 0. In the ordered phase 
with r < T c the leading contribution to the longitudinal 
correlations is given by [66j | 



(-iri(s-)i 5 



(24) 



where \(S*)\ is the staggered magnetization correspond- 
ing to the onset of long range antiferromagnetic order as 
given by Eq. (fT5j) . In contrast, in the disordered phase 
with r > T c , the correlations decay exponentially with a 
power law prefactor (r36l | 



(-iy 



-1/4 



(25) 



where Z" 1 = ln(r/r c ). To confirm this behavior we de- 
scend deep in to the Mott phase in order to suppress 
the effects of the small XY terms given by Eq. (jl~4")) . As 
shown in Fig. [10] the results arc in excellent agreement 
with the theoretical predictions ([2"4"1) and (f2"5j) . 



VIII. ATOM-MOLECULE CORRELATIONS 

The above considerations are consistent with the no- 
tion that away from the Ising transition line the con- 
nected correlations decay exponentially in the MI. As 
advocated in Ref. [32j the absence of low lying gapless 
excitations precludes the possibility of the novel supcr- 
Mott behavior proposed in Ref. f39l [4p . This gained 
further support in a recent comment 4l| which confirms 
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the presence of exponential decay in the atom-molecule 
correlation functions. We discuss these observations for 
the bosonic Hamiltonian ([3]) within the framework of the 
Ising description. 

In order to investigate the possibility of countcrflow su- 
percurrents [(34|, [65| of atoms and molecules, the authors 
of Rcf. Hdl consider the connected correlation function 



0.02 



C am (n) = (TO t (n)a(n)a(n)a t (0)a t (0)™(0)) 

- (m t (n)a(n)a(n))(a t (0)a t (0)m(0)). 



(26) 



Using the spin mapping ([5]) [32l . |33| this may be written 
in the equivalent form 

C am {n) = 2 [<S+(n)S-(0)> - (5+(n)) (S~ (0))] . (27) 

Deep within the second Mott lobe one may therefore use 
the Ising Hamiltonian @ to gain a handle on the bosonic 
correlation function (|26|) . To this end we may decompose 
the first term in Eq. (j2"7| as 

(Si~S i+n ) = (SfSf +n ) + (SfSf +n ), (28) 

where the mixed component terms cancel. It follows that 

C am (n) = 2 [{S?S? +n ) - (S*){S? +n ) + (SfSf +r7 )} , (29) 

where we use the fact that (Sf) = for the transverse 
field Ising model (Q| in zero magnetic field. That is to 
say, C om (n) is the sum of the connected xx correlation 
function and the yy correlation function. From our previ- 
ous discussion in Sec. IVIIl these contributions generically 
decay exponentially, in agreement with the findings of 
Ref. [4l|. However, along the locus of the Ising quan- 
tum phase transition one expects power law contribu- 
tions. Using Eqs. ([T8"]) and (fT9")) one obtains 



(30) 



C am (n) ~ 2 [(27m)- 2 - (-1)" (B / '4)n~ E '< /4 



In order to confirm this characteristic oscillatory depen- 
dence we perform DMRG on the ID bosonic Hamiltonian 
([3j with L = 96 and periodic boundaries. Employing the 
finite size result (|20[) obtained by conformal transforma- 
tion one obtains 



Cam {n) 



1 



tt/L 



sin I 



- (-1)" 



B 



tt/L 



sin I 



9/4 



(31) 



As shown in Fig. I 111 the results are in excellent agreement 
with the underlying Ising correlation functions. 



IX. CANONICAL SOFTCORE BOSONS 

Having provided a discussion of the model ([3]) in the 
reduced Hilbcrt space, with r a = 2 and r m — 1, we turn 
our attention to the more general problem with canoni- 
cal softcore bosons. In this situation one must allow for 
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FIG. 11. Correlation function C am (n)/2 at and away from 
criticality obtained by DMRG on the ID system Q with 
L — 48 and periodic boundaries. We use the same parameters 
as in Fig. [5] (a) At criticality clear oscillations are present in 
conformity with the Ising description. The solid line corre- 
sponds to the theoretical prediction in Eq. (|31l) . (b) Away 
from the critical point with V = 2r c , C am {n) shows a finite 
correlation length due to the Ising gap. 



the presence of additional virtual intermediate states in 
the magnetic description. For example, in a configura- 
tion with two pairs of a-atoms on neighboring sites the 
softcore problem allows virtual hopping to take place, 
in contrast to the problem with r a = 2. This modifies 
the coefficients, J zz , h, and C, but the Ising description 
remains valid deep within the second Mott lobe. One 
again obtains the effective magnetic Hamiltonian given 
in Eq. Q but with the modified coefficients, 



Jz 



4f 



Un.m, U n 



4f 



Un 



and 



2e a - U aa + 



mi 

U aa 



At, 

Urn' 



(32) 



(33) 



together with 



C = L — 



Uaa 

2 



-it, 



r 



8 Un 



U a 



Uar 



mi 

Uaa 



At 2 



U„ 



(34) 



as shown in Appendix [X] Note in particular that the ef- 
fective longitudinal magnetic field, h, now depends on the 
hopping parameters and is therefore already present at 
second order in the strong coupling expansion. In order 
to see the effect of these additional intermediate states it 
is instructive to examine the change in the ground state 
energy of the bosonic Hamiltonian ([3]) upon increasing 
the local atomic Hilbcrt space restriction, r a . For sim- 
plicity, we consider U mm oo, corresponding to hard- 
core molecules with r m = 1. As shown in Fig. IT2T a). 
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FIG. 12. Ground state energy density of the Hamiltonian Q 
obtained by DMRG with L = 20 (to aid comparison with 
ED results for the Ising model) and periodic boundaries, (a) 
We use the parameters in Fig. [5] and increase the local atomic 
Hilbert space restriction r a , with r m = 1 held fixed. The lines 
are results for the energy density of the Ising Hamiltonian l|4"|) 
obtained by ED (L = 20) with J zz = 5.02 x 10~ 4 and h = 
for r a = 2, and J zz = 4.23 x 10 -4 and h = 7.89 x 10 -5 
for r a — 3, 4. The change from r a — 2 to r a = 3 is due 
to the presence of additional virtual states which modify the 
Ising model coefficients. The absence of any change beyond 
r a — 3 is consistent with second order perturbation theory 
around the second Mott lobe, (b) Ground state energy density 
for the same parameters as in panel (a) with the additional 
interaction U m m = 4. We increase the local Hilbert space 
from r a = 3 and r m = 2 to r a = 4 and r m = 3. The absence of 
any further change is consistent with the maximum occupancy 
for virtual states explored in second order perturbation theory 
around the second Mott lobe. The solid line is the energy 
density of the Ising Hamiltonian @ obtained by ED (L — 20) 
with J zz = 4.16 x 10" 4 and h = 7.27 x 10" 5 as given by Eqs. 
(|32l) and (|33[) respectively. 



the ground state energy changes on going from r a = 2 
to r a = 3 due to these additional hopping processes. 
However, increasing the atomic restriction beyond r a = 3 
has no further effect, since higher occupations are not ex- 
plored at second order in perturbation theory within the 
second Mott lobe. The excellent agreement of the bosonic 
results with r a = 3 and r a = 4 therefore supports the 
applicability of the second order Ising description. The 
solid lines shown in Fig. [l~27a) correspond to the ground 
state energy density of the Ising model ((4]) obtained by 
exact diagonalization with the appropriate coefficients. 
Our DMRG results are in excellent agreement with these 
predictions. As shown in Fig. H"27 b) this agreement also 
extends to the generic softcore problem. At second order 
in perturbation theory around the second Mott lobe the 
maximum local occupation is r a = 3 and r m = 2. The 
absence of any further change in the ground state en- 
ergy on increasing the local Hilbert space justifies both 
the use of this truncation and the Ising description. In 
Fig. [T3] we show the existence of an antiferromagnetic 
Ising quantum phase transition within the second Mott 
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FIG. 13. Entanglement entropy difference AS(L) = 
Sl(L/2) — Sl/2{L/4) showing an Ising quantum phase tran- 
sition within the second Mott lobe for restricted softcore 
bosons. We truncate the local Hilbert space to r a = 3 and 
r m = 2 which is the maximum occupancy explored at second 
order in perturbation theory. We use the same parameters as 
in Fig. I12f b) and set t = 0.005. The peak height corresponds 
to a central charge c w 1/2. 



lobe for these restricted softcore bosons. The entangle- 
ment entropy difference as given by Eq. (fl"2|) exhibits a 
peak on transiting through the magnetic transition. At 
criticality AS = | In 2 + . . . and the peak height is in 
good agreement with an Ising quantum phase transition 
with central charge c = 1/2. 



X. FERROMAGNETIC ISING TRANSITION 

An interesting aspect of the canonical softcore result 
(|32p is that one may explore both antiferromagnetic and 
ferromagnetic Ising interactions due to the relative mi- 
nus signs. This is readily seen by exact diagonaliza- 
tion on small systems using the techniques employed in 
Refs. [32|, [33|, |63|. As shown in Fig. [T3]an Ising transi- 
tion indeed persists with ferromagnetic parameters and 
h = 0. This is also confirmed by DMRG results for the 
ground state energy of the bosonic Hamiltonian with 
L = 64 and periodic boundary conditions. In the ther- 
modynamic limit, the ground state energy density of the 
Ising model ^ is given by [66|, Eq/L = C/L + eao, where 

1 

^00 — ; 

47T 



dky/4T 2 + Jl z + 4TJ ZZ cos k. (35) 



Defining E' = Eq — G the result for Eq/L is in good 
agreement with Eq. ([33]) as shown in Fig. [15] (a). In 
addition the finite size corrections to the ground state 
energy are well described by the conformal result [zO, [7l| 



7TCV 

6X2 



+ 



(36) 



where v = \J zz \/2 is the effective velocity and c is the 
central charge. Going beyond our h = parameter 
choice, this opens up the exciting possibility of study- 
ing the celebrated Es mass spectrum of the ferromag- 
netic Ising model in a magnetic field (46| . As exempli- 
fied by recent experiments on the quasi-one-dimensional 
Ising ferromagnet CoNF^Og (cobalt niobate) [U EU, it 
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FIG. 15. (a) Rescaled ground state energy E' = Eq — C of 
the bosonic Hamiltonian ((3| with the FM parameters used 
in Fig. 1141 The DMRG data are in good agreement with 
Eq. (|35|) for the thermodynamic limit of the Ising model, (b) 
The residual finite size corrections, E' /L\J ZZ \ = eoo/|J zz | — 
7tc/12L 2 + . . . , are consistent with the central charge c w 
0.499 of an Ising transition to a FM state in the absence of a 
longitudinal magnetic field. 



FIG. 14. ED results for the Hamiltonian ([3]) with hardcore 
molecules (r m = 1) and up to three atoms (r a = 3) per site. 
We set e a = 0, U aa = 2, U am = 4, t a = 2t m and take t a = 0.01 
corresponding to FM exchange with J zz ~ —3.94 x 10 -4 . For 
simplicity we set h — e m — 2 + Qt\ = by taking e m = 1.9994. 
(a) The rescaled energy gap Ai = E\ — Eq shows a clear 
intersection at T c ~ 1.969 x 10 -4 ~ 0.4997|J ZZ | corresponding 
to a FM transition, (b) Scaling collapse with the Ising critical 
exponent, v — 1. (c) The rescaled pseudo-magnetization m — 
(| Sf\)/L indicates a transition at the same value of T c . 
(d) Scaling collapse with the Ising critical exponent /3 — 1/8. 

would be interesting to probe the dynamical correlation 
functions of the ID bosonic Hamiltonian ([1]) in more de- 
tail. The non-trivial excitations will manifest themselves 
at specific frequencies dictated by the emergent Eg mass 
spectrum. Similar behavior is also expected in the AFM 
case in the presence of a staggered longitudinal magnetic 
field. 



XI. CONCLUSIONS 

We have investigated the Mott insulating state of 
bosonic pairing Hamiltonians using analytical and nu- 
merical techniques. We have described the behavior of a 
broad range of physical observables, including local ex- 
pectation values and correlation functions, within the 
framework of the paradigmatic quantum Ising model. 



As advocated in Refs. [lH, HH, Hl[ our results are con- 
sistent with the absence of super-Mott behavior within 
the second Mott lobe. Our results differ from the usual 
two-band Bose-Hubbard model which exhibits counter- 
flow superfluidity [6J, |65| since Feshbach resonant pairing 
interactions favor a distinct Mott phase with either two 
atoms or one molecule per site. As such, XY terms gener- 
ically appear at cubic order in the hopping strengths and 
are intrinsically suppressed. An alternative way to see 
this is to note that for finite Feshbach coupling, the sym- 
metry of the Hamiltonian is reduced from U(l) x U(l) 
down to U(l) x Z2. As such, one naively expects Ising 
transitions to occur in a Mott phase with large phase fluc- 
tuations and an unbroken U(l) symmetry. Nonetheless, 
one cannot rule out the possibility of novel transitions in 
other regions of the phase diagram due to higher order 
terms in the strong coupling expansion becoming appre- 
ciable. It would be interesting to explore this in more 
detail. In addition, we highlight the possibility of using 
these systems to explore the E 8 mass spectrum of the FM 
Ising model in a magnetic field. There are many direc- 
tions for further research including the influence of higher 
bands and the possibility of realizing Ising transitions in 
Bose-Fermi mixtures; see Appendix |B| 
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Appendix A: Effective Magnetic Hamiltonian 

To obtain the effective spin description within the sec- 
ond Mott lobe we perform degenerate perturbation the- 
ory. We partition the Hamiltonian as H — Hq + V where 



H a = J2 WK 1 - Pot) + ( e ™ ~ h a /2)V Ql 



(Al) 



and 



v = -J2 tvQLfii* + H - c -) + h f+ h o s * Voi - ( A2 ) 

(ij)a i 



Here h = e m - 2e a - U aa , V Ql = \$) . . + and 

(A3) 



ttI _ \ ^ , \ ~* Uaa' 



This approach is appropriate deep in the Mott lobe and 
in the vicinity of the transition where g, ho and t 2 /U are 
small and V may be treated pcrturbativcly. Up to second 
order in V the effective Hamiltonian is 

H cS = V (h + V + V{1 - Vo) 1 - (1 - V )V ) V 

(A4) 

where Eq = L(e m — Hq/2) and Vo — Yli^oi- Since the 
hopping terms are the only source of coupling between 
the degenerate subspace spanned by Vo and the remain- 
ing Hilbert space, the second order term reduces to 

H$> = £ t 2 a V blb ja 1 b] a b ta V + H.c. (A5) 



{ij)a 



This acts on two sites so within the degenerate subspace 

Vo, we must consider its action on four basis states. The 
(2) 

action of H^ s on sites with neighboring molecules yields 

4£ 



H 



Uram + ^0 ' J 

and for an arrangement of neighboring atoms gives 

12t? 



(A6) 



H. 



(2) 



'eff W/i W/j - TT , \V/i l-VV j 

When the neighboring species arc different 



(A7) 



2t;, 



ft 

' ■■■>■ 



U a rn U a a 2 U a Ti 



(A8) 



In order to obtain the second order contribution to the 
effective Hamiltonian we collate these terms. In addition 
we exploit the spin projection identities Vo = 

(l/2 + S?)V and (^Vo = {l/2-S*)P and expand 
the resulting expression to leading (zeroth) order in h : 



H™ ~ Vo J zz J2 S * S j + h *Y, S * + C 2 ] Po, (A9) 
where J zz is given by Eq. (|32[) and 



z (Ylti 4^ 



2 \ U aa Ujnr, 



(A10) 



The constant offset is given by 



zL 



C 2 = - — 



8 \U, 



4f 



12i 



4i 



U aa U a 



(All) 

Inclusion of the zeroth and first order contributions in V 
yields the effective Ising Hamiltonian §4§ with h = ho+hz 
and C = Eo + C 2 as given by Eqs. ([33]) and (J34J). In 
order to obtain the effective Hamiltonian with r a = 2 
one must exclude the contribution from Eq. (|A7[) as the 
intermediate states involved in the perturbation do not 
satisfy the imposed constraint. Similarly if r m = 1 the 
contribution from Eq. (|A6j) must be excluded. 



Appendix B: Bose— Fermi Hamiltonian 

Throughout this manuscript we have focused exclu- 
sively on the bosonic homonuclcar Hamiltonian (jTJ) and 
the associated Ising description. However, it is evident 
from the general setup shown in Fig. [T] that similar re- 
sults may also emerge with more than one atomic species. 
For example, this is confirmed in Refs. [HI, for the 
heteronuclear bosonic case. In this appendix we note 
that a similar Ising description may also apply with two- 
component fermionic atoms and bosonic molecules. We 
consider the Bosc-Fcrmi Hamiltonian 

H = ^2e a n ia -^2^2t a (d\ a d ja + H.c.) + H F 



ia (ij) a 

X 1 Uaa' x U mrn . . 

i,a^£a' i 



(Bl) 



where a =4-,t are fermionic atoms, a = m is a bosonic 
molecule, n ia = d\ a d la and H F = gJ2i( d l m d il d n + ]1 - c -)- 
Assuming the existence of a second Mott lobe as depicted 
in Fig. [TJ with either two fermionic atoms or a bosonic 
molecule per site, one again obtains an Ising Hamiltonian 
© acting on the states |J|) = d|d||0) and \ = d^O). 
The parameters in Eq. (j4]) are now given by T = 2g and 



Jz 



t 



t 



2/ 2 



2 Uf m + U± m Ul m — Ufl Uf m — Uf± U mm 

(B2) 
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The effective magnetic field is given by 



h = e m -(e l + e t + U u )-—^, (B3) 

U rnm 



and the constant offset is 

°- £ M + t-£)- (B4) 

In the limit U mm -> oo the results (|B2t . ([B3]) and ([B4]) 
coincide with the results of Sec. IIIII where we identify 
E% = u aa and C/j. m = L^m = Vara- In view of our 
findings in the bosonic problems it would be interesting 
to examine the Bose-Fermi mixture (|J31|) in more detail. 
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